Scalar and vector form factors of the in-medium nucleon 
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Abstract 

Using the quark-meson coupling model, we calculate the form factors at a- and w-nucleon strong- 
interaction vertices in nuclear matter. The Peierls-Yoccoz projection technique is used to take 
account of center of mass and recoil corrections. We also apply the Lorentz contraction to the 
internal quark wave function. The form factors are reduced by the nuclear medium relative to 
those in vacuum. At normal nuclear matter density and = 1 GeV^, the reduction rate in 
the scalar form factor is about 15%, which is almost identical to that in the vector one. We 
parameterize the ratios of the form factors in symmetric nuclear matter to those in vacuum as a 
function of nuclear density and momentum transfer. 
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The change of hadron properties in a nuclear medium is of fundamental importance in 
understanding the implication of QCD for nuclearphysics. One of the most famous nuclear 
medium effects may be the nuclear EMC effect and it has stimulated theoretical and 
experimental efforts to seek nuclear quark-gluon effects for almost two decades. 

Recently, the search for modification of the electromagnetic form factors of bound protons 
has been performed in polarized (e, e'p) scattering experiments on ^^O and ^He nuclei j2|. 
The experiments measured the ratio of transverse to longitudinal polarization of the ejected 
proton, which is proportional to the ratio of electric to magnetic form factors of a proton. 
However, conventional calculations including free-proton form factors, appropriate optical 
potentials and bound-state wave functions as well as relativistic corrections, meson-exchange 
currents (MEC), isobar contributions and final-state interactions, fail to reproduce the ob- 
served results in ^He Indeed, full agreement with the experimental data was only 
obtained when, in addition to the standard nuclear calculation^a change in the form factors 

Recent inclusive neutrino experiments on ^^C at Los Alamos ^5] also suggest that the mea- 
sured total cross section is about a half of the standard, relativistic shell model calculation 
including final-state interactions within the distorted wave impulse approximation 6]. In 
the neutrino reaction, the charged- current vector form factors of bound nucleons are slightly 
enhanced, while the axial form factors are quenched by the nuclear medium 0]. Finally, the 
effect of the bound nucleon form factors reduces the total cross section by about 8% relative 
to tKat ca.c„.ated witK tKe .ee fonn facto, g We stress that t>.s oo„.eot,o. is caused by 
the change of the internal quark wave function at the mean-field level and hence there is 
no obvious double counting with MEC etc. This is a new effect which should be taken into 
account additionally to the standard nuclear corrections. 

Furthermore, the measurements of polarization transfer observables in exclusive (p, 2p) 

n n 

proton knockout reactions from various nuclei [21, 110| again indicate that it is difficult to ac- 
count for the measured polarization transfers within the conventional, relativistic distorted 
wave impulse approximation [l^]. To reproduce the measured spin observables, it is nec- 
essary to simultaneously reduce the scalar (a) and vector {uj) coupling constants and the 
meson masses by about 10 ~ 20% In particular, the analyzing power {Ay), polarization 
(P) and spin transfer coefficient {Dgs') are very sensitive to the change of a- and tu- nucleon 
coupling constants and their masses in a nuclear medium. These may again imply the change 



2 



in the internal structure of bound nucleons. 

If the quark substructure of the nucleon is modified depending on the nuclear environ- 
ment, it would leave traces in a variety of processes and observables, including various form 
factors. These modifications of bound nucleons can be successfully described within the 
context of the quark-meson coupling (QMC) model jl^. In the model, the medium effects 
arise through the self-consistent coupling of a and uo mesons to confined quarks, rather than 
to the nucleons. As a result, the internal structure of the bound nucleon is modified by the 
surrounding nuclear medium. 



The electromagnetic form factors o 



proved cloudy bag model (ICBM) 
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bound nucleons jj] have been studied using an im- 



IJ, together with the QMC model. In the ICBM, a 



simplified Peierls-Thouless projection technique (the weight function w{'p) appearing in the 
nucleon wave function is assumed to be unity) is used to account for center of mass (cm.) 
and recoil corrections. In addition to it, a Lorentz contraction of the internal quark wave 
function is included. The axial form factor in nuclear matter has also been calculated in a 
similar manner Furthermore, the form factors at a- and tu-nucleon strong- interaction 
vertices in a nuclear medium should also be investigated. The change of these form factors is 
very significant in understanding how the strong interaction is modified in nuclear matter. It 
is also expected to play an important role in analyzing the polarization transfer observables 
in the exclusive (p, 2p) reactions js , 10 1 . 

In this Letter, we study the scalar and vector form factors at a- and tu-nucleon strong- 
interaction vertices in symmetric nuclear matter. We shall calculate these form factors using 
a relativistic constituent quark model with a harmonic oscillator (HO) la] or a linearly rising 
(LR) confining potential 16| and the Peierls-Yoccoz (PY) projection technique. If we use the 
"minimax" principle (or the saddle point variational principle) it is easy to obtain an 

approximate solution to the Dirac equation with any potential. Since we choose a Gaussian 
wave function for a confined quark as ansatz, it is possible to calculate the form factors 
analytically and thus transparent to see how the PY projection and the Lorentz contraction 
of the quark wave function work in the form factors. Instead, in this exploratory study, we 
do not include the pion cloud effect which can explicitly be treated in the ICBM. (We will 
study this effect in a forthcoming paper.) 

In the QMC model, the mean-field approximation is applied to the a and uj meson fields, 
which couple to confined {u or d) quarks in nuclear matter. Each quark then satisfies the 
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Dirac equation 

[-la ■ V + + t/eo„/(r)]^(f) = E,iP{f), (1) 

where m* = — g^a and Eg = tg — g^u with the quark energy. We take the free quark 
mass rrig to be 300 MeV. The mean-field values of a and u mesons are respectively denoted 
by a and uj, and g^ and g^ are the corresponding quark and meson coupling constants. We 
use a confining potential of HO type, Uconfi^) = (c/2)(l + /57°)r^, or a LR one, Uconfi^) = 
(A/2)(l + l3'~f^)r, where (3 {0 < [3 < 1) controls the strength of the Lorentz vector- type 
potential. The potential strength is taken to be c = 0.04 GeV^ or A = 0.2 GeV^ jisj ]. 

Although for the LR potential the Dirac equation cannot be solved analytically, the min- 
imax principle allows us to obtain an approximate solution very easily and accurately jl7|. 
Since the Dirac Hamiltonian does not have a lower bound for the energy spectrum, the usual 
variational method cannot be applied. The minimax principle amounts to minimizing (max- 
imizing) the energy expectation value of the upper (lower) component of the quark wave 
function with respect to variational parameters. A trial wave function for the lowest-energy 
state is usually chosen as 

^(^1 = A^of.,,51. (2) 
\it,ba ■ ru[r) ) 

with A^o a normalization constant, M(r) = g^^^*"^/^ and h and ^ the variational parameters. 
These parameters are determined so as to minimize the quark energy tg with respect to h 
and maximize it with respect to ^. Note that for the HO potential with /3 = 1 this gives the 
exact solution jigj ]. 

First, we fix the parameters of the model in vacuum. The nucleon mass in vacuum 
(a = = 0) is given by Mtv = 3eg — eo, where eo accounts for corrections of cm. and gluon 
fiuctuations. The parameter eo is fitted so as to obtain the free nucleon mass M^[= 939 
MeV). The minimax principle then determines the parameters h and ^. These values are 
given in Table HI 

In matter, the scalar field couples to the confined quark and hence the quark mass changes 
depending on the nuclear environment. The nucleon mass in matter is then reduced 
because the a exchange induces an attractive force between nucleons. In an iso-symmetric 



nuclear matter, the total energ}[_ 



expression in the QMC model 



per nucleon) at nuclear density is given by the usual 
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Etot = r + ^-a^ + ^PB. (3) 

pBK^T^r J ^PB 2m^ 



where mo-(= 550 MeV) and m^(= 783 MeV) are respectively the a and u) meson masses, 
and gu}{= 3(7^) is the cj-nucleon couphng constant. The values of the scalar and vector 
mean fields are, respectively, determined by self-consistency conditions: {dEtot/da) = and 
{dEfot/duj) = 0. The latter condition ensures baryon number conservation, while the former 
gives a transcendental equation for the scalar field in matter. 

The coupling constants are fitted so as to reproduce the nuclear saturation property 
[Etot — Mn = —15.7 MeV) at normal nuclear matter density po(= 0.17 fm~^). Note that for 
each value of one has to use the minimax principle to obtain the in- medium parameters 
b and ^. The coupling constants and nuclear properties at po are listed in Table HI The 
cr-nucleon coupling constant g^^ is defined in terms of the quark scalar density Sjy: = 
Sg^Sjy^d- = 0), where Siy{a) = J dr %l){r)ip{r). 

The wave function for a nucleon moving with momentum p can be constructed by the 



PY projection technique 
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^(n, ra, fa; p) = N{p} J dx e'^'^^n, fg, fg; f), (4) 
where N{p) is a momentum- dependent normalization constant 

[Nip)r'= /dfe-^>lp(f)]^ (5) 



with 

Pirl = (6) 

Here (f) is the quark wave function in momentum space. The localized state $ is simply given 
by a product of the three individual quark wave function 

<l>(fi, r2, fa; x) = ■?/'(fi - x)ip{f2 - x)^p{fs - f), (7) 

where x refers to the location of the center of the nucleon and fj {j = 1, 2, 3) specifies the 
position of the j-th quark. 

Because the nucleon consists of three point-like quarks, the expectation value of an opera- 



tor with respect to the nuc 
quark expectation values 



eon wave function Eq. (j3)) may be given by a sum of the individual 
. In the Breit frame, where the initial (final) momentum of 
the nucleon is taken to be —q/2 {q/2) with q the momentum transfer, the scalar and vector 
form factors are respectively given by 

r(.)(g2) = ^[N{Q')f j dfe^'l■^^i,{T)i^\i,{T,q), (8) 



where = —q^ + = q'^, and we ignore a small tensor term at the tu-nucleon coupling. 
ip in Eq.(jH]) is represented by 

dk 



m <1)=J J^s e''-'mW{k, q), (9) 
where 

W{k, q) = Jdf e-'(^^/2+^)"lp(f)]^ (10) 



and 



P(r) = l^s e-''-'mm. (11) 
Now we can calculate the scalar and vector form factors in nuclear matter analytically: 

r(.) (Q , Ps) - e ^^^yz^ . (12) 



where = Q'^/b'^ 



and 



z,(e) = 1 + + + ^e, (13) 
Me) = j^jf . (16) 
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n^-^ia = i-i^J^'+Uk'^ii^'' ^^^^ 



1 ; 32 ^ 128 ^ 1152^ ' ^ ^ 

= ±^e. (19) 

Recall that the variational parameters ^ and b (thus x^), which appear in the quark wave 
function, depend on ps- We have renormalized the vector form factor so that r„ = 1 is 
maintained at zero momentum transfer. The scalar form factor is also rescaled by the same 
factor as in the vector case 

In contrast, if the cm. correction is ignored, the form factors are given by 



1 + 



(20) 



Because is small and ^ ^ 0.5 (for Pb/pg < 2.0) at small momentum transfer, we can 
expand the form factors. Up to 0{x^) or 0{^'^), we find that r°^^ = 1 — {^^^'^ — while 
Eq. lfT^ gives F^^^^ = 1 — Q^,^^ — /Q. The cm. correction thus moderates the reduction of 
the form factors. 

Apart from the cm. correction, it is also vital to include the Lorentz contraction of the 



2l|. The full 



internal quark wave function at moderate or large momentum transfer 
form factors F^^s-j can be obtained through a simple rescaling 

F..^(Q2) ^ V,s.{Q') = r^T..^(r/*Q2), (21) 



where rj* = {M^/E'^Y with E^f = yMj^ + <5^/4. The scaling factor in the argument arises 
from the coordinate transformation of the struck quark and the prefactor r/* comes from the 
reduction of the integral measure of two spectator quarks in the Breit frame Thus, 
the scaling factor rj* (in vacuum t] with Mn) should appear in any nucleon(baryon)-meson 
form factors if the nucleon (baryon) is assumed to have a three-quark cluster structure. 

To illustrate the effects of the cm. correction and Lorentz contraction on the form factors, 
we show in Fig. ^ the vector form factor in vacuum. The cm. correction considerably 
enhances the form factor in comparison with the result without both effects (see the dotted 
and dot-dashed curves in the figure). The effect of Lorentz contraction is also important. 
If the Lorentz contraction is ignored, the form factor drops away like ~ ^-^'^/^ at large 
Q^. The inclusion of the Lorentz contraction removes this objectionable exponential falloff. 
Because of the factor t], the form factor is proportional to 1/(1 + Q^/A^) and is modified 
to x^/(l + (5^/A^) with A = 2Mn (see also Eq. pTj) ). As a result, the inclusion of the Lorentz 
contraction enhances the form factor at large (see the dot-dashed and solid curves). 

Because our aim is to study the density dependence of the form factors in nuclear matter, 
we consider the ratios of the in-medium form factors to those in vacuum: 

f..^(Q^p5) 

R{^\ {Q , Pb) = ■ (22) 

^ F^.)(Q2,pB = 0) ^ ^ 

The form factors in symmetric nuclear matter F^s^ are thus given by 

F(s.{Q\pb) = R(s^{Q\pb) X FfT.'m, (23) 



where Ff^^ are the form factors empirically determined in vacuum 2M ■ -^ig- El the ratio 
of the in-medium scalar (vector) form factor to that in vacuum is illustrated as a function 



of and ps- (Because the ratios for the LR potential are similar to those for the HO 
potential, we focus on the HO case for a while.) At ps/po = 1 and = 1.0 GeV^, the 
in-medium scalar (vector) form factor is reduced by 15 (14)% relative to that in vacuum. 
The reduction rate depends on (3 very weakly. By contrast, at Pb/Po = 2 and = 1.0 
GeV^, the scalar form factor decreases by 35 (29) [24]% for /5 = (0.5) [1.0], while the 
vector form factor diminishes by 28 (26) [22]% for /? = (0.5) [1.0]. At high density the 
dependence of the reduction on (3 is thus rather strong, and the reduction rate is correlated 
with (see Table HI). 

As in the case of vacuum (see Fig. Q), the effect of Lorentz contraction is again seen at 
large Q^. For example, at Pb/Po = 2 and = 1 GeV^, the vector form factor with the 
Lorentz contraction is about 7% larger than that without it. We also note that, in the HO 
case with j3 = 0.5, the full vector form factor gives the root-mean-square radius of 0.53 fm. 
If we neglect the Lorentz contraction effect, it is 0.46 fm. 

Finally, we parameterize the ratios for the scalar and vector form factors in Eq. p3|l . 
Such parameterizations are very useful in analyzing the experimental results, e.g., for the 
exclusive (p, 2p) proton knockout reactions With an error less than 0.2%, the ratios 

can be represented by 

R(^s-^{Q\pb) = I + A^s-^iQ^)ipB/po) + B^.^{Q')ipB/po)\ (24) 



where 



^0.06829\ /0.2302\ , /0.1845\ ^ /0.04613\ 3 

^ " V0.04072/^ ' 

^^^y^ = 0.3738)^+ (0.2668 - (0.06494)^ ' (2^) 



^^^^^ " "Uo6323j " U2464j^+ Ul71lj^' " ^04072^^'' ^^^^ 
'0.3856\ , /0.3021\ 2 /0.07763\ 3 



and 



0.005071\ /0.02499\ /0.09473\ 2 /0.1070\ 3 /0.03914\ 



^'^^^ l^0.003569j 1^0.03304/^^1^0.1167;^ 1^0.1245^^^1^0.04474/^'^^'^^ 
/0.04296\ /0.2081\ 2 /0.2380\ 3 /0.08967\ . , , 

= -Uo550oj^+ U227lj^ - (0.2509;^ + Uo9337j^ ' ^^'^ 

In Eqs. fl25|) ^ ()28|) . the upper (lower) numbers are for the case of the HO (LR) potential with 
P = 0.5 (1.0), which provides the effective nucleon mass M^/Mjy = 0.71 ~ 0.72 at po (see 
Table nj. The in-medium form factors are thus given by Eqs. (j23j) and (j24j) . 
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In summary, using the QMC model we have calculated the form factors at a- and u- 
nucleon strong-interaction vertices in symmetric nuclear matter. We have applied both the 
PY projection technique and the Lorentz contraction of the internal quark wave function. 
The form factors are reduced by the nuclear medium relative to those in vacuum. The 
cm. correction moderates the reduction of the form factors in matter, and the Lorentz 
contraction is vital at large momentum transfer. We have found that the reduction in 
the scalar form factor is about 15% at Pb/po = 1 and = 1 GeV^. This rate is almost 
identical to that for the vector form factor. In contrast, the scalar and vector form factors are 
respectively reduced by about 30% and 25% at p^/po = 2 and = 1 GeV^. The reduction 
of the form factors is expected better to reproduce the polarization transfer observables 
measured at RCNP and iThemba laboratory 0, We have parameterized the ratios of 
the form factors in symmetric nuclear matter to those in vacuum. This provides a convenient 
formula to estimate the in-medium form factors. It is very intriguing to re-analyze the data 
of polarization transfer observables for exclusive {p, 2p) proton knockout reactions 
including the modification of both the form factors and meson masses in matter . 
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TABLE I: Coupling constants, eo, b, ^, and nuclear incompressibility K. The parameters eo, 
b and ^ arc fixed in vacuum, wliilc b* , ^* and are calculated at normal nuclear matter density. 
Here cq, b and K are quoted in GeV. The value of P is specified in the parenthesis in the first 
column. 
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FIG. 1: Vector form factor in vacuum (for the HO case with /5 = 0.5). The fuh result is denoted by 
the solid curve, while the dot-dashed curve shows the result with the cm. correction but without 
the Lorentz contraction. The result without both corrections fEa. (|2U() ) is denoted by the dotted 
curve. 





FIG. 2: Ratios for the scalar (left panel) and vector (right panel) form factors in the case of the 
HO potential with /? = 0.5. 
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